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KUilSRIGAL EVALUATIOlf OP ^L'HE c-IKTSGRAL OCCURRUK} HI 
TEE THS0DGR3SH ARBITRAiTr AIRFOIL POTENTIAL TIEIORY 

By Irven Haiman 



A more precise method of evaluating the c-integral 
occurring in the arbitrar-r airfoil theory of Theodorsen 
(IfACA Reps. Fos. 411 aiid 452} is developed by retaining 
hJ.gher order temc i;i the Taylor expansion and by iise 
of Siiapaon's rule. Forrsulas are ^iven for routine calcu- 
lation of the e -Integral and foi* the necessary coapu- 
tctiont-l coefficients. The computational coefficients 
are tabulated for a 40-polnt division of the range of 
integration fron 0 to Sir. ?:ith no increase in conpu- 
tational work the syatomatic error in the nunerical 
value of e is reduced from the order of 1 percent to 
approximately 0.1 percent. 



nJTRO!>UCTIOH 



The solution of the general problem by means of 
confonnal transformation for the florr about an arbitrary 
airfoil (references 1 and 2) , a. ?7?nietrical lattice 
(reference 5), and a biplane (reference 4) involves the 
determination of the inaginai^y part of a conplex trans- 
formation function, given the real part. As shor^n in 
references 1 and 2 the real part may be expanded in a 
Fourier series and the imaginary part is the conjugate 
Fourier series. It is also shovm in these references 
that the imaginary part e may be obtained from the 
real part \!/ by the following functional equation; 



e(cp«) = / >!'(<P) cot ci<P (1) 



2 
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Tbls Integral occurs frequently in confomal- 
transforraation problens Involving the evaluation of the 
functions on the circle. A procedure for the wanerical 
evaluation of this integral is given in references 1 
and 2. This method, trhich is currently in use at IIIAL, 
gives an error of about 1.5 percent for a 40-point 
division of the range of lntegre.tlou. An improvement 
in the accuracy is therefore very desirable, particularly 
if the labor involved is not increased. A revised method 
fiven herein is found to involve a little less work than 
that previously req\iired and to give an error of only 
about 0.1 percent. Constants for lis e in this more pre- 
cise method have been coaputed and are presented In 
table I. 



EyALTIATICH OP TH3 C-I3JTI3GRAL 



The evaluation-' of the e-integral is complicated by 
the discontinuity at cp = (pf . This difficulty nay be 
surmounted by a separate solution across the discon- 
tinuity* \7hen s = cp - (p' is substituted in eqi;ation (1) 



T r2Tr-(?« 

€(<?') = -j^ / *(<P' + s) cot f ds 



or, because of the periodicity of this function. 



:(<p«) =-^f / + s) cot y 



ds (2) 



The discontinuity noT' cccirs at s = 0. For purposes of 
numerical evaluation chis Intejrcl 3iay be broken up as 
f ollous : 

ifr^ r2if-B "I 

£(cpi) o -;iJ I i|r(<p> + a) oot |.ds + / t{Q>' + a) cot ids 
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Where 



«1 = 



r3 

-Fir I 
•J- a 



+ s) cot ^ ds 



(4) 



and 



'2 



Jr2Tr-a 
s 



« + s) cot ^ ds 



(5) 



Evaluation of €^.- The first liitegral includes 

the dlscontlnxilty and the limit s s.ay be taken as some 
convenient small value. By a ^fi^v-lor series expansion the 
Integral Is easily evaluated as follows: 



\|/(<p» + s) = \K<?') + s^^'(qp') + -^i|/"((p«) + |^\1/'"(<P') + 



V.hen this expansion is substituted in equation (4) the 
integrals containing the even-ordered derivatives are 
found to be identically zero. Equation (4) then becomes 



1 



J-s 



a cot ^ ds + 



3 S 

s cot ^ ds + 



rhere the derivatives are evaluated at <p' . The Taylor 
expansion for cot ^ is 



a. 3 _ 2 s 3 

cot ^ - - - ^ - ^ - ... 
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and equation (6) Is thus obtained as 



4. 



Evaluation of ^g.- Hhe second Integral eg (equa- 
tion (E)) nay be rearr&iiged for convenience In numerical 
calculation as follows: 

= -^jj^>K<5>' + s) cot I ds+ jf - s) cot I dsj 



= -^fP&C®' + a) - iir(«?t - a)] 



cot ^ ds 



(8) 



where -s has been substituted fcr s in tb© pccond 
Integral and the limits have beer, reatrangecl accordingly. 



f . .r.-.V^ 5 - ^l?!?-.! • " I'^J^ T-cxci-sroe 1 tt^ fpterval 

2Tr/n (n is? fr. o-'^n rtvnrcp'5r-) □ v -?Vaf5 oi" . re 

desiCiV'.tfsd v V p ..i. ^-j. ^V. 

where "^q Is the value of ^ at (p = (p» ana >!/ = 
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is the -value at <P = <P' ± Tf. The Integrations are per- 
formed over intervals of width 2v/n with the >1/ values 
at the midpoint of the interval. The range of integra- 
tion for en la froia s = -ir/n to s = ir/n and for 
ep from ir/n to Sir- -2,. 

The first integral e^^ is evaluated by retaining 
only the first-order terms in s^ 

€i = -gVs^o' = -1*0' (9) 



where the slope ^ITq' is determined graphically at >|fQ, 
that is, at <p = 9' (s = 0). 

The second integral is composed of the sum of 

integrals across each interval 



2k+l 




\|f(<p« + s) cot I ds 



The function \U does not change nuch across the 
Interval and is therefore approximated by use of the 

value at the midpoint s t|(^(p» + —^j. Then, 



n-l .^^^ 
^2 --^^^K 



k=l 



n Q 

cot ^ ds 



2k-l 

-=1T 
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or, "bj equatlcn (8), with ^^jj^ = ^^j^ 



n/g 
lc=l 



(10) 



Td^ere 



ajj. = log 



Sn 



fin 5 0 



(11) 




(12) 



Values of the constants ajr were given in reference 1 

for n = 10 and in reference 2 for n = 20. Revised 
values for these constants, together rrith those for 
n - 40, are given in table II. 

Improved nethod .- The numerical accuracy of the 
evaluation of the € -Integral will be shomi to be 
larproved by the follovring method: The interval 0 to 2v 
is divided into n equal parts and the ^ values are 
designated as in the previous section. The second 
integral €3 is evaliiated by Simpson's rule from 

to ■>l'n«i (^n-1 ~ ^-l)' range of integration for €] 

is therefore tvrlce as large as that in the previous sec- 
tion, that is, s -Str/n to 2Tr/no The approximation 
in which only the first-order term of equation (7) is 
used is insufficient and the higher derivatives must be 
used. These derivatives are most conveniently obtained 
by numerical differentiation. 
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The ITe'^tcn-Stirllnc forE:ulasfor derivatives (refer- 
ence 5y p. 75) are 



(15) 



where s = tabular interval (2Tr/n) . ?he piean central dif- 
ferences 6* can te expressed in tcrLns of the tabular 
values as 



2SM/ = 


^1 








2G^>lf = 


*2 


- Sxlr^^ + 


2^-1 - >l'-2 


>(14) 


Z&H = 


*3 


- 4*2 + 


5*1 - 5*_i + 4*_2 - *_3 





gives 



The substitution of relations (13) in equation (7) 



+ S6 



5.713 s' 



450 1512 



- )- 1 

• •Of"* aaer 



(15) 
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The further substitution of relations (14) in equa- 
tion (15) gives 



or 



«1 - ^l(*-l - ^l) + ^2(^-2 - *2) + ^3(^-3 + *3) tlS) 



vrhere 



„ 1 /239 167 „2 • \ 

- Msir - ^52^ ^ ~ '"J 

^2 = -gir [(I Ms) " ^^(2^ sfe) - • • •] 
-1/88 2 ^2 



> (17) 
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The second Integral €g Ig evaluated by Sinpson's 
rule from ilf-j^ to il'n-i 

'2 = -^i[*l <^ot g ^ 4.^2 cot f 



+ 2*5 CDt ^ + . 



. + cot 



or, by equation (8), v;it;i = >!r_ic 



n/2 



(13) 



where 



01 5 ^ H 



ci: = 35 coo 



cn = 0 



(19) 



where (e^ccept in the first term) o = 2 for k odd and 
a = 4 for k even. The conplete inte^^ral is given hy 
« = «1 + ^2* or 




(20) 
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where 

A]j = + Cjj (21) 



Values of Ajr for n = 40 are given in table I. 



ACCURACY OP FVALUATIOH 



The accuracy of the two methods of evaluation 
described nay be determined by integrating variou'i 
harmonlcp. The results are presented as ratios of the 
integrated value to thd correct value eo that a value 
of unity is a correct evaluation. Values of this ratio 
for the haraonlcs ere: 



Harmonic 


40-polnt method 
of reference 1 


40-point present 
method 


1 


1.01434 


0.99977 


2 


1.02662 


1.C0044 


3 


1.03700 


.99940 


4 


1.04547 


1.00060 


5 


1.05280 


. 99340 


10 


1.06112 

■ 


. 95420 



Inasmuch as the higher harmonics generi:lly enter 
in a much smaller proportion than the lower ones for 
such contours as are encountered for airfoil shapes, the 
error of the 40-polnt method of reference -1 is of the 
order of 1.5 percent, nvhereas that of the nen method 
presented herein is approximately 0.1 percent. 
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TABLE I.- VALUES OP A^ PGR USE WITH EQUATIOIT (20) 











K 




k 




n = 40 


i 


n = 40 


1 


0.52827 


t 

11 


0.01423 


2 


.14824 


12 


.02422 


5 


.07614 


13 


.01C21 


4 


.10259 


14 


.01698 


5 


,04024 


r 1? 


.00690 


6 


.0G542 


1 16 


.01083 


7 


.02720 


17 


.00400 


8 


.04588 


18 


.00528 


9 


.01S51 


19 


.00066 


10 


.0?533 


20 


0 



TABLS II.- VALUES 0? aj^ PCS USE NITH ECJUATIOIT (12); 
IIETEOD OF E3TKaSTO33 1 AlID 2 











k 












• 


n = 10 


n = 20 


n = 40 


1 


1.06544 


■ 1,09057 ' 


1.09666 


2 


.44311 


•4942e 


.50671 ■ 


3 


.23117 


.31141 


.33028 ■< 


4.,- 


.10302 


,21750 '•• 


.2^303 ^ 


5 


0 


r .1577-5 


.19028 ■ 


6 




.11440 ■■ 


. 15453 


7 




■ c0802-: 


.12341 


8 




. ,05115 . 


.10026 


9 




■' ,02493 1 


.09207 ? 


10 




0 


.C78G2 I'- 


11 






.06714 " 


ir; 






.05711 


15 






.04 316 - 


14 






.C4004 


15 






.03255 


IC 






.02553 


17 






.01337 ■' 


18 






.01245 ' 


19 






.OColG f 


SO 






0 



